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The nonlinear interaction of DC current flowing in a thin metal film with an external low-frequency 
AC electromagnetic field is studied theoretically. The nonlinearity is related to the influence of 
the magnetic field of the DC current and the magnetic field of the wave on the form of electron 
trajectories. This magnetodynamic mechanism of nonlinearity is the most typical for pure metals 
at low temperatures. We find that such interaction causes sharp kinks in the temporal dependence 
of the AC electric field of the wave on surface of the sample. The phenomenon of amplification of 
the electromagnetic signal on the metal surface is predicted. We also calculate the nonlinear surface 
impedance and show that it turns out to be imaginary-valued and its modulus decreases drastically 
with the increase of the wave amplitude. 



I. INTRODUCTION 

As well-known metals possess quite peculiar nonlinear electrodynamic properties (see., for example, Refs.EUl). In- 
deed, nonlinearity in a response of plasma or semiconductors to an electromagnetic excitation is usually associated to 
significant deviation of the electron subsystem from equilibrium. To the contrary, nonequilibrium in metals as a rule 
is weak due to large concentration of charge carriers. Nevertheless, a nonlinear regime in these media is rather easy 
to observe. Such a situation is possible because of the fact that nonlinearity in metals does not necessarily related 
to the nonequilibrium of the electron subsystem. Nonequilibrium is caused by the existence of a weak electric field, 
while nonlinearity is originates from strong magnetic fields. The Lorentz force, determined by the magnetic field of an 
electromagnetic wave or magnetic field of the transport current, affects the dynamics of charge carriers. As a result, 
the conductivity of a sample placed in an AC electromagnetic field depends on the spatial distribution of magnetic 
field of the wave. Such a magnetodynamic mechanism of nonlinearity is inherent to pure metals at low temperatures, 
if the mean free path of conducting electrons is rather large. 

Magnetodynamic nonlinearity causes a number of nontrivjal phenomena in the electrodynamics of metals. As an 
example, one can mention generation of the current statesau in a sample placed in a DC external magnetic field. 
The sample acquires a DC magnetic moment if irradiated by an additional strong AC electromagnetic field. The 
magnitude of the magnetic moment depends in a hysteretic manner on-an external DC magnetic field. Under current 
states conditions, a hysteresis-like interaction of radiowaves is observecH as well as the appearance of electromagnetic 
dissipative structuresfj. This specific mechanism of nonlinearity in metals causes a decrease of collisionless damping of 
heliconso. Therefore, the spiral w.aves with large amplitudes can x>ropagate even in conditions_when there is no their 
linear electromagnetic excitationsQ. Magnetoplasmic shock wavesu and soliton-like excitationsE2l are also predicted for 
the regime of strong magnetodynamic nonlinearity. 

In the present paper, we study a novel manifestation of magnetodynamic nonlinearity, namely, interaction of an 
external electromagnetic wave and a strong DC transport current in a thin metal film, which is also displayed in a 
quite unusual way. The sample of thickness d is assumed to be much smaller than electron mean free path I, 

d<Z, (1) 

and electron scattering on a surface of the film is supposed to be diffuse. It is known0 that in the static case (when 
external AC signal is absent), the magnetic field of a current can essentially affect the conductivity of a thin metal 
specimen and, thus, its current-voltage characteristics (CVC). In this situation the value / of the current is rather 
small so that the typical radius of curvature R{I) of electron trajectories in magnetic field is much greater than the 
film thickness, 

d<£R(I), R(I) = cp F /eH(I) cx J" 1 . (2) 

Here — e and pf represent electron charge and Fermi momentum, respectively. In Ref.0, it was shown that nonlinear 
peculiarities of CVC are connected to the antisymmetrical spatial distribution of the magnetic field of the DC current 
inside the sample. The magnetic field equals to zero at the middle of the film and takes the values H to —H at the 
opposite boundaries of the sample, where 



H = 2-kI/cD. 



(3) 
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In this formula, c denotes the speed of light in vacuum and D is the sample width. Spatially alternating field of the 
DC current entraps a part of electrons in a potential well. Trajectories of such particles are flat curves winding around 
the plane of alternation of the magnetic field. The relative part of the trapped electrons in the order of magnitude is 
equal to the typical angle (d/R) 1 ^ 2 <C 1 of their crossing of that plane. Taking into account that trapped carriers do 
not collide with the film boundaries and interact with the electric field along their whole free path I, one can write 
the following estimating formula for their conductivity atr'- 

a tr ~ a^d/Ril)) 1 ' 2 oc I 1 ' 2 . (4) 

Here <jq represents the conductivity of the bulk sample. At the same time, there exist flying electrons which do 
collide with the boundaries of the specimen and, according to Ref.__, have the conductivity of the order of ao(d/l). 
Apparently in the range of rather strong currents, when the inequality 

(dR(/)) 1/2 « I, (5) 

holds, the conductivity of the film is determined by the group of the trapped carriers. As a result, we observe the 
deviation from the Ohm's law: the voltage U is proportional to the square root of current, 

U oc I 1 ' 2 . (6) 

For the film with thickness d = 10 -3 cm, the electron mean free path I = 10 _1 cm, and Fermi momentum pp = 10 -19 
g- cm/s the nonlmearity_becom.es noticeable ((dR) 1 / 2 ~ I) at values of the magnetic field H(I) about 1 CLe. The 
theory developed in Ref.EJ is in a good qualitative agreement with experimental data (see, for example, Ref£_). 

In an external magnetic field h, collinear to the magnetic field of the current, the plane of the sign alternation for the 
magnetic field shifts to one of the two boundaries of the film (see Fig. 1). That in turn leads to appreciable diminution 
of the number of the trapped particles and, therefore, their conductivity. In particular, such a situation would take 
place under symmetrical irradiation of the film by the low- frequency electromagnetic wave of large amplitude. The 
frequency is supposed to be so small that the AC magnetic field h(t) of the wave is virtually uniform across the metal 
(i.e. the wave penetration depth 8 is much greater than the sample thickness d). Then the conductivity of metal 
essentially depends on time and, therefore, strong nonlinear effects in the sample response to the AC electromagnetic 
excitation should appear. Being of interest from the both theoretical and experimental points of view this problem 
has not been investigated yet. 

In the present paper we study theoretically the temporal dependence of electric field at the surface of the film, 
which carries a strong DC current of the fixed value I satisfying inequalities (g) and (g). It is shown that with an 
increase of the amplitude h m of the AC magnetic field this dependence becomes anharmonic, turning into a series 
of sharp nonanalytic peaks. The case of sufficiently high amplitudes h m > H, when the total magnetic field in the 
sample is spatially alternating during some part of the wave period and has constant sign during the other part, is 
of particular interest. In such a situation, the electric field is also characterized by kinks in its temporal dependence 
due to the periodical appearance and disappearance of the group of trapped carriers. The effect of amplification of 
the electric signal on the film surface is predicted as well. It turns out that, because of the presence of the strong 
DC transport current in the sample, the absolute value of the AC electric field of the wave is //(di?) 1 / 2 ^> 1 times as 
many as the corresponding magnitude in the absence of the DC current. 

We also calculate the nonlinear surface impedance of the film, which turns out to be pure imaginary value in the 
main approximation in the parameter d/5 <C 1, and show that its modulus monotonically diminishes with the growth 
of the AC amplitude decreasing in If(dR) 1 / 2 ^> 1 times. Simultaneously the conductivity of the trapped particles 
falls down and, consequently, S increases. 

II. PROBLEM STATEMENT AND GEOMETRY 

Consider a metal film of the thickness d with a DC current I flowing along. The sample is irradiated from both 
sides by a monochromatic electromagnetic wave with a magnetic component collinear to the magnetic field of the 
current. The x-axis is oriented perpendicularly to the film boundaries. The plane x = corresponds to the middle of 
the sample (see Fig. 1). The y-axis is directed along the current, and the z-axis is parallel to the vector H(x, t) of the 
total magnetic field which is a sum of the DC magnetic field of the current H(x, t) and the AC magnetic field of the 
wave h(x, t), 



H(x, t) = {0, 0, H(x, t) + h(x, t)}. 



(7) 
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FIG. 1: Geometry of problem. Schematic representation of the trajectories of trapped (1), flying (2), and surface (3) electrons. 



The film length L (the dimension along the y-axis) and its width D (z-axis dimension) are much greater than the 
sample thickness d. We assume diffuse scattering of the electrons on the him boundaries. Maxwell's equations in the 
assumed geometry can be written as 



dH(x,t) _ 4tt dE(x,t) l dH(x,t) 

dx c dx c dt 



(8) 



where j(x,t) and E(x, t) represent the y components of the current density and the electric held. Boundary conditions 
for Eqs. (N) are 



H(±d/2, t) = h m cos ut =F H. 



(9) 



Let H be the absolute value of the magnetic held on the surface of the metal him and h m denote a wave amplitude. 
According to Eq. (||), the held H is determined by the DC current /. No special relation between magnitudes H and 
h m is assumed. 

We consider a quasistatic situation when the wave frequency to is much less than the relaxation frequency v of the 
charge carriers, 



u) <C v. 



(10) 



Here we suppose that the AC magnetic held inside the sample is quasiuniform and virtually does not differ from its 
value on the sample surface, h(x, t) ~ h rn cos cut. In other words, the typical spatial scale S(ll>) of variation of the AC 
magnetic held is much greater than the him thickness d. Furthermore we assume that the curvature radius R(x, t) of 
electron trajectories in the total magnetic held Ti.(x. t) is also much greater than d, 



rf<<5(w), d<^R{x,t), R(x,t) = cp F /e\H(x,t)\. 



(11) 



III. ELECTRON DYNAMICS, CURRENT DENSITY, AND CVC OF FILM 



Let us consider electron dynamics in the nonuniform AC magnetic held H.(x,t). We shall assume the following 
gauge of the vector potential: 



A(x, t) = {0, A(x, t), 0}, A(x, t)= dx'H{x', t) 



(12) 
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FIG. 2: Phase space (p y ,x). Regions of existence of flying (I), trapped (II), and surface (III) particles in spatially sign alternating 
(a) and of constant sign (b) total magnetic field. 

It is suitable to choose the lower limit of integration in Eq. (|l^) depending on whether or not there exists the plane 
x = Xo(t) of the sign alternation of the magnetic field TL{x,t) at the present moment. This plane exists during the 
time intervals when h m | cos u>t\ < H because the values h m cos ujt — H and h m cos cut + H of the total magnetic field at 
the film boundaries have opposite signs (see Eq. (0)). In this case, one should take xo(t) as the lower limit in integral 
(|l2|). Then the vector potential A(x,t) is negative. It reaches its maximum value (which equals to zero) at the point 
x = xo(t). Within other time intervals, when the inequality h m \ cosu>t\ > H holds, the magnetic field Ti.(x,t) inside 
the sample is of a constant sign. In such a situation, one should choose sign (cos ut)d/2 (sign(x) is the sign function) 
as a lower limit of the integration. In this case, vector potential, also being negative, vanishes at one of the boundaries 
of the film. 

The integrals of motion of electron in the field H(x,t) are the total energy (it equals to the Fermi energy) and 
the canonical momenta p z = mv z and p y — mv y — eA(x,t)/c (m is the electron mass). The electron trajectory in a 
perpendicular to the direction of the magnetic field plane is determined by the velocities v x (x,t) and v y (x,t). In the 
case of Fermi-sphere with radius pf = mv, we obtain 

\v x {x,t)\ = (vl -«jj) 1/a , v± = {v 2 -v 2 z ) 1/2 , v y (x,t) = (p y + eA(x,t)/c)/m. (13) 

Classically allowable regions of the electron motion along the x axis are determined by the inequalities, 

— p y — mv± < eA(x,t)/c < —p y +mv±. (14) 

These inequalities provide the positivity of the radicand in Eq. (|l3| ) for \v x (x,t)\. 

The regions of the electron motion in the phase plane {x,p y ) are described schematically in Fig. 2 for two cases: 
when there exists the plane x — xo(t) of the sign alternation of the magnetic field 7i(x,t) (Fig. 2, a) and when such 
plane is absent (Fig. 2, b). For definiteness we have chosen the moment of time when the magnetic field of the wave 
is positive (costot > 0). The upper border on the phase plane is described by the curve p y = mvj_ — eA(x, t)/c and 
the lower one is given by p y = —mv± — eA(x, t)/c. The electrons are naturally divided in three groups depending on 
the sign and value of the integral of motion p y . Below, we give inequalities determining the regions of their existence 
at an arbitrary moment of time. 

1. Flying electrons: 

Py = — ttwx — eA[— sign(cosa;£)d/2, t]/c < p y < mv±, \x\ < d/2. (15) 

These particles collide with the both boundaries of the film. Their trajectories do not twist significantly because 
of d <C R(x, t). Flying electrons exist at every moment of time irrespective of the presence of the plane x — xo(t) 
(i.e irrespective of the relation between h m cos ujt and H). 
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2. Trapped electrons: They appear during the periods of time when h m \ costal < H and the total magnetic field 
H.(x,t) within the sample passes trough zero. Their states are bounded by the region (see. Fig. 2, a), 

—mvx < P y < Py = —mv± — eA[sign(cosujt)d/2, t]/c, 

a;*(£)sign(cosa;i) < icsign(cosa;i) < d/2. (16) 

Here a;* (t) represents the breakpoint of the trapped electron most distant from the film boundary. One can find 
it from the equation, 

A{x t ,t) = A[sign(coscjt)<2/2,*]. (17) 

According to Eq. (|l6|), this electron group occupies the region x*(t) < x < d/2 when cos u>t > and the region 
—d/2 < x < x*(t) if cos u>t < 0. The trajectories of trapped particles are almost flat oscillating curves due 
to periodical motion of the particles along x-direction and the uniform motion along the y-and z-axes. The 
temporal period of oscillations with respect to the plane x = xq equals to 2T, where 

dx 

T= / , , .. . (18) 

J xi (t) MM) 

The breakpoints X\(t) and x%{t) (x\(t) < Xo(t) < X2(t)) are the roots of the equation, 

eA(x li2 ,t)/c= -mv± - p y . (19) 

3. Surface electrons: 

These particles collide only with one of the boundaries of the film. In our case of diffuse scattering of the 
electrons on the surface, their influence on the nonlinear conductivity of metal is neglig iblefcl Thus, we do not 
take them into account thereafter. 

The current density of the flying and trapped particles can be deduced by means of solving the Boltzmann kinetic 
equation. One should linearize the kinetic equation with respect to the electric field E(x, t), which can be represented 
as a sum, 



E(x,t) = E + £(x,t), 

(20) 



1 f dA(x,t) dA(t)\ 

£M = - c {—t dT) 



Here the first term, Eq, is a potential (uniform) component and £{x, t) is a rotational (nonuniform) field of the wave. 
Spatial averaging of the latter over the x-axis direction gives zero. The value A(t) represents a spatially averaged 
magnitude of the vector potential, 

1 rd/2 

A(t) = - / A(x',t)dx'. (21) 
dJ-d/2 

The magnetodynamic nonlinearity is accounted for in the kinetic equation by means of terms which contain the 
total magnetic field TC(x, t) = H(x,t) + h(x, t) entering the Lorentz force. We calculate the current density in the main 
approximation with respect to the small parameter d/6(u>) (see Eq. (0)). In this approximation, as it was mentioned 
above, the AC magnetic field h(x, t) becomes spatially uniform and is equal to its boundary value, h(x,t) = h m costvt. 
The electric field is also independent of the ^-coordinate and coincides with thp-.value Eo(t). For the case of uniform 
electric and external magnetic fields, the current density was obtained in RefO. If the conditions (|^) and (||) hold 
the following asymptotics for the current density of the flying and trapped electrons are valid: 

jfl(t) = CTfl(t)E (t), 

a f l & = l ao J ln d^ , R ±( f ) = c PF/e\h m \ cos ujt\±H\, (22) 



jtr{x,t) = <T t r{x,t)E (t), 
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367T 1 / 2 re i i / z 

a tr{x,t) = 5r U) a °\c — [ A ( x,t ^ ~ A (s [ gn-{cosu}t)d/2,t)} J , 



1/2 



5r 2 (i/4) 



■CPF 



(23) 



x»(i)sign(cos wi) < a;sign(cos cut) < d/2. 

In the limit u> — > 0, Eqs. (22) and ( p3| ) transform into the corresponding formulae of Ref.0. 

Let us substitute the current density in the first of Maxwell's equations (||) for its asymptotic expressions 
( p3| ) and introduce a dimensionless coordinate and vector potential, 

£ = 2xsign(cos u>t)/d, a,(£,i) — A(x, t) /A(sign(cos cut) d/ 2, t). 

The equation for the quantity a(£, t) has the form, 

r[l-a(£,f)] 1/2 + l, &(*)<£<!, 
1, -!<£<£*(*), 



9£ 2 



|) and 

(24) 

(25) 



£*(t) = 2x sr {t)sign(cosujt)ld. 

The dimensionless coordinate confines the region of existence of the trapped particles and, accordin 
and (p4]), satisfies the equation, <z(£*, £) = 1. The parameter r represents the ratio of the maximum maj 
conductivity of the trapped electrons to the conductivity of the flying particles, 



(26) 



a tr (x ) 967T 1 / 2 / 



\A(sign(cosuit)d/2, t) 

CPF 



1/2 



u n 5r 2 (l/4)d 

The dimensionless quantity, u, is related to the voltage U = EqL on the sample 

U 



]n~\R+/d). 



g to Eqs. |H 
nitude of the 



(27) 



cL|A(sign(coswt)ri/2, t)\/wafid 2 
Equation ( p5j ) should be solved together with the boundary conditions, 

da(l,t) d h m \cosujt\—H 



(28) 



2 A(sign(cosuJt)d/2,t) , 



da(-l,t) 



I COSLUtl 



H 



2 A(sign(cos wt)d/2,t) ' 



a(l,t) = 1. 



(29) 



and the third one is a consequence of 



The first two of these expressions are dimensionless boundary conditions 
normalization ( p4]) of the vector potential. 

Within the interval < £ < 1, the solution of Eq. (^5|) is symmetrical with respect to the point £o(t) = 

(1 + £*(i))/2, where the dimensionless vector potential reaches its minimum value (which equals to zero, a(^Q,t) = 
<9a(£o, £)/<9£ = 0). This solution is described by the formula, 



|e - £ (*)| = (3/4™) 1 / 2 r i6,t d([l - (1 - C) 3/2 + 3C/2r]- 1 / 2 
Jo 



(30) 



One can not obtain the field distribution and the current density within the region of existence of the trapped electrons 
in an explicit form. However, by means of Eq. ( BO]), i t is possible to calculate the average magnitude of the conductivity 
of the trapped carriers (|23|) within the intervalp™, 



<Jtr_ 



r / dC(l-C) 1/2 [l-(l-Cr /2 +3C/2r]- 
Jo 

(^C[l-(W) 3/2 + 3C/2r]- 1/2 ) ■ 



-1/2 



(31) 
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The bar above a tr denotes spatial averaging. In the remaining region of the sample (—1 < £ < £*(i)), there exist only 
flying electrons, and the solution of Eq. (]2^) is given by the following formula: 

a(& t) = 1 - (2u) 1 / 2 (l + 2r/3) 1 /2( ? _ + „(£ _ ^)) 2 /2. (32) 

Expressions (|3(i|) and ( |32"| ) and their derivatives are sewn together at the point £ = The solution given by 

Eqs. ( |30| ) and (|32|) contains three parameters, £o: u an d r, which should be found from boundary conditions ((2^). It 
is essential that the value A(sign(coscj£)<i/2, t) of the vector potential appearing in Eq. ( p9| ) is not an independent 
parameter due to its relation to r via formula (p7|). 

Adding term by term the first two boundary conditions in Eq. ( p9f ) and using Eqs. ( p0[ ) and (32), and (|28|), we find 
the following expression for the drift of the plane x — xq: 

£o = 2a;osign(coswt)/d = Lh m \ cos Ldt \ ^ ^ cosujt\ < H. (33) 

2irU<jfid 

In order to determine the value of u (i.e. the voltage U), let us integrate the left and right-hand sides of Eq. fl25| ) 
from -1 to 1 taking into account the boundary conditions for the derivative da(t;,t) / d£ in (E9h. The integral of the 
function [1 — a(£, t)] 1 ^ 2 appearing in the right-hand side can be reduced to the product 2(1 — £o)&tr/rcrfi with the use 
of the condition a(l, t) = 1. Taking this into consideration as well as formulae ( pS| ) and ( p3| ) for the quantities u and 
£o, we have after some simple transformations, 

,, cL H(I) + (a tr / a f i)h m \ cos ujt\ , , ., . „ 

2ird<Tfi{t) l + atr/a-fi 

According to Eq. (^lf), the ratio of conductivities, atr/cfi, depends on the parameter r. Using expression (g£J) for it, 
relation ( p7[ ) between A(sig~D.(cosu>t)d/2, t) and r, and solution (|30|), we obtain from the first boundary condition in 
Eq. ( p^ ) the algebraic equation for r, 

i? — /i m | cosu>t\ \ 2 C/ 



r(1+2r/3) H — i — ■) m^jiy (») 

Here we have introduced the following notations: 

^ = 25r 4 (5/4) c PF rf - _ AclLH 
9n el 2 ' 37rcrorf 2 

The parameters H and [7 represent those magnitudes of the magnetic field and voltage for which the characteristic 
length (Rd) 1 / 2 of the arch of electron's trajectory is of the order of the mean free path I. 

Expressions (|3l]), ([54|), and (^) define, in an implicit form, dependence of the voltage U on the current I for the 
case h m \ cosuit\ < H . At these conditions there exists the plane of the alternation of sign of the total magnetic field 
within the sample. If the opposite inequality, h m \cosujt\ > H, is valid, the trapped electrons are absent (r = 0, 
£* = 1, <j tr = 0) and CVC is described by the formula, 

U = 2 ^ ( ^ , h m \costut\>H. (37) 

As seeing from formula (|34|), the voltage on the sample displays nonanalytical behavior vs. time: the dependence 
U(t) has kinks at the moments when the AC magnetic field h m cos u>t vanishes. This is an essentially nonlinear effect 
caused by the contribution of a large group of trapped electrons into the electric current. 

The temporal dependence of voltage ([34j) for the case when the wave amplitude is not too large (h m < H) and there 
exist the trapped carriers during the whole period 2tt/lu is shown in Fig. 3, a. Fig. 3. b represents the dependence U(t) 
for the opposite case h m > H, in which during some part of the wave period (at h m \ coswtj > H) the conductivity is 
caused by the flying particles only. 



IV. NONANALYTICAL TEMPORAL DEPENDENCE OF ELECTRIC FIELD 



Knowing the vector potential A(x,t), one can calculate the rotational electric field E(x,t) as a correction to E (t), 
(see Ref. (p0|)). We are interested in the difference AS (t) — £(d/2,t) — £(—d/2,t). This value is proportional to the 
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FIG. 3: Time dependence of the voltage U at relatively small (a, h m < H) and large (b, hm > H) AC amplitudes. 



rate of alteration of the magnetic flux trough the cross-sectional plane, which is perpendicular to the direction of the 
vector of the total field 7i(x,t), and thus can be measured in experiment. 

From Eqs. (|3(i| ) and (|32"1), it follows that the difference a(l,i) — a(— l,t) is connected to the derivatives da(l,t)/d£ 
and da(—l,t)/d£ by the relations, 



o(l,t)-a(-l,t) = -&(t) 
a(l,t) - a(-l,t) 



da(l,t) da{-l,t) 
da(l,t) da(-l,t) 



h m \ cos(jjt\ < ff, 
/i™ I cos ujt\ > H. 



(38) 
(39) 



Let us now turn to the dimensional variables in Eqs. ( p8[ ) and (39) using boundary conditions (|29j) and relation 
( p7j ) between the values of A(sign(coswt)<i/2, t) and r. After that one can obtain the following expression for the 
magnitudes of the vector potential at the film boundaries: 



A(sign(cosujt)d/2,t) = - H d\n {R+ / d)r 2 /4, 



A(-sign(coswi)d/2,t) = - H d\n {R + / d)r 2 / A ~ 2H\x (t)\ 



(40) 



at 



hm,\ cos ujt\ < H 




and 



A(sign(cosu;t)d/2,t) = 0, A(— sign(cos ut)d/2, t) = — dh m \ coswt] 



(41) 



at 



h m \ cos u>t\ > H. 



Formulae ( |40| ) and ( |4l| ) are sewn at the time moment when h m \ cosmt] = H. The parameter r in Eq. ( |27| ) vanishes, 
and the plane x — xa(t) coincides with one of the boundaries of the sample, |xo(£)| = d/2. From relations ( [40| ) and 
(pp|), by means of formula j33|) for £o(i), we derive the expression for the difference A£(t) of magnitudes of the electric 
field at the film boundaries, 



A£(t) 



2Hdx (t) H(I)Lh m d 



Of 



2tt 



dt 



COS Ojt 



a fl (t)U(t) 



h m I COS LOt\ < H. 



(42) 



If the inequality h m < H holds, the previous relation is valid during the whole period of the wave. However, in the 
case h m > H, there exists a time interval when the plane x = xo(t) of alternation of the sign of the total magnetic 
field is absent. If such a situation takes place one should use formula (|df) in order to obtain the dependence A£(t). 
Finally we come to the result below, 



A£(t) = A£l sinwt, A£l = dh m uj/c, h m \ cosujt\ > H. 



(43) 



From this, it follows that the difference A£(t) is a harmonic function of time, i.e. the response of the film on the 
external electromagnetic excitation turns out to be linear if there are no trapped electrons. It is obvious that formula 
([l3|) also describes the dependence A£(t) at small magnitudes of the current / (H <C H), when the contribution of 
trapped particles to the conductivity is negligible during the whole period of the wave. Then the value A£l represents 
the amplitude of a linear response. 

The dependence A£(t) is shown in Fig. 4 for a wide range of the AC amplitudes h m and for the large magnitudes 
of the DC magnetic field H of the current /, when the inequality H ^> H (or inequality (||)) is valid. It is obvious 
that the ratio of the amplitude A£ m to its linear value A£l does not depend on h m . From relations ([42|), (|34|), and 
(E35|) at cos cvt = 0, we find the expression for A£ m , 
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The ratio AS m /ASL is determined by the magnitude of the DC magnetic field H and can be much greater than 
unity. In other words, there exists an effect of amplification of the electric signal at the film surface. For small AC 
amplitudes (curve 1, h m — H/300) the signal turns out to be quasi-harmonic. However, with the increase of h m 
the dependence AS (t) shows kinks. Curve 2 has kinks at the points of extremum, i.e at the time moments when 
the AC magnetic field h m cosojt vanishes. These singularities are related to the nonanalytical behavior of CVC of 
the film (see. Eq. j34|) and Fig. 3). Curve 3 corresponds to the case h m — 5H/3, in which the trapped electrons 
are absent during a part of the wave period. In such a situation, the dependence AS (t) contains additional kinks 
arising at the moments of appearance and disappearance of the plane x = xo(t) of the sign alternation of the total 
magnetic field. They are located symmetrically with respect to the points of extremum as shown in curve 3. By 
means of formulae p^), d43|) , ([34)) , and (|35|), we find the right and left derivatives of the function AS(t) at the point 
t = (1/uj) arccos(i? / 'h m ) of the first kink, 



d AS{t) 



dt AS L 



t=t -o 



ujH 

hm. 



(45) 



d AS(t) 



at as l 



t=t +o 



ujH 

hm 



2\n(R+/d) \H 



hV' 2 fh 2 



H 2 



(46) 



According to Eq. (46), the right derivative is negative and has large absolute value even at [(h m /H) — 1] > 1. 



SURFACE IMPEDANCE OF FILM 



Let us analyze the dependence of the surface impedance at the film boundary x — d/2 on the AC amplitude h m 
under conditions of interaction of the transport current and the electromagnetic wave. The impedance is proportional 
to the ratio of the first Fourier harmonics of the electric S^ and magnetic h^, fields at the surface of the sample, 



47T Su, 



8nS^ 
c h„ 

2ttc 



„ _ fdA(d/2,t) OA 



2-kc 



dt 



dt 



e iuJt dt 



2-k/u 



(A(d/2,t) ~A(t)) e iujt dt. 



(47) 



Taking into account Eqs. ([53]) and ([l0]), we deduce the boundary value of the vector potential for the periods of time 
given by the inequality h m \ coswtj < H, 



A(d/2,t) = 



-Hd\n 2 (R + /d)r 2 /4, atcostut > 0, 



-H\n 2 (R + /d)r 2 /4 : + cHLh m cosujt/2irU(t)<jfi{t), atcosojt < 0. 
In the case h m \ cosu;t\ > H, the following expression is valid (see. Eq. (f4l|)): 



A(d/2,t) = { 



0, 

dh rn COS Ujt, 



atcosujt > 0, 
atcosust < 0. 



(48) 



(49) 



Let us calculate the mean value of the vector potential A(t) for h m \cosLot\ < H, when there exists the plane of 
alternation of sign of the field. According to Eqs. (pG) and (p2|), we have 



A(t) 



A(sign(cos u;t)d/2, t) 



-(2/3)u(t)e 



a&tn = &(*) + (2«(t)) 1/2 (l +r(t)/3) l ' 2 i 2 {t) 
-l 

1/2 



Ar{t)u(t) 



(d( 



o y/1 - (1 - C) 3/2 + 3C/2r(i) 



(50) 



In the case h m \cosujt\ > H, one should use solution (|32|) with r = 0, ^* = 1 in order to find A(t). Proceeding to 
dimensional variables and using Eqs. (^8|) and (|37|) , one can easily obtain 



A(t) 



dh rn I COS LUt\ 



-Hd. 



(51) 
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We draw reader's attention to the fact that the mean value of the vector potential depends on time only via the 
term |cosa;£|, A(t) — A(\cosuit\). This follows from formulae (|35|), (p8|), and (33) for the values r , u , and £ as 
well as from the relation (27) between A(sign(cosu;t)d/2, t) and r. It also implies that the surface impedance in the 
main approximation with respect to d/S has imaginary part (reactance) only. The latter is a consequence of the full 
transparency of the film. 

We start calculation of the reactance X with the case of relatively small amplitudes h rn < H, when the group of 
trapped electrons exists during the whole period of the wave. Let us substitute expressions ( |48| ) and ( |50| ) into Eq. p7|). 
Then, the integrals containing A(t) and —Hd In 2 (R + /d)r 2 /4 vanish since these functions depend on | coswi| only. By 
means of formula (|34|) for the voltage U, the remaining integral can be transformed into the form, 



X 



8duj 

~~~2~ 



n/2 



l + a tr {T)/cr f i{T) 



1 + {<J tr {r)/ a fi{T)){h m /H) cos-. 



cos 2 rdr, hm < H. 



For the case of large amplitudes h m > H, one should calculate the reactance using formulae 
It represents a sum of two terms, 



(52) 

J), and f3). 



X 



8dw 



r/2 



1 + a tr (r)/afi 



7r/2-arcsmff/h m 1 + Kr (t)/<7/j (t) ) {h m /H ) ( 



■ cos 2 rdr 



+ 



7r/2— arcsin H/h n 



cos rdr 



at h m > H. 



(53) 



The first term corresponds to the temporal interval when the trapped electrons exist in the sample, and the second 
one is related to the interval when these particles are absent. 

Let us calculate the asymptotics of the surface reactance for the case of rather large amplitudes h m ^S> H. For this 
purpose, we rewrite integral (p3|) in another form, 



X = 



8duj 



ir/2 



cos rdr 



o 



zr/2 

7r/2— arcsin H/h^ 



1 + <Jtr{T)/(Jfl 



1 + {a tr {r)/afi(T))(h m /H) cost 



1 ) cos 2 rdr 



In the second integral, we substitute the variable of integration (h m cos t)/H 
power series in the ratio H/h m . Then one finds, 



(54) 



r\ and expand the integrand in a 



X 

xl 



i + -(H/h m y 

IT 



l + <j tr {ir/2)/a fl (ir/2) 
l + a tr {ir/2)/<T f i(Tr/2)r 1 



r) drj, 



where 



2vr 

X L = —ud 

c z 



(55) 



(56) 



is the same as the value of reactance in the absence of the DC transport current. The conductivities at r (^/2), and 
<t/z(7t/2) are taken at the moment of time when the AC magnetic field h m cosujt turns into zero. Therefore, their 
ratio is much greater than unity due to inequality (Q. Taking into account condition (]44|), we calculate integral ( p35| ) 
and obtain the following asymptotics for the reactance, 



X 

X~2 



2 

3^ 



H 



H < h„ 



(57) 



Now we consider the case of the extremely small amplitudes described by the inequality h m <C 
H<jfi(ir/2)/<7 tr (Tr/2) ~ (HH) 1 / 2 . The integrand in Eq. ( |52] ) can be presented as a power series in /i m /(i?ff) 1 / 2 . 
As a result the asymptotic takes the form, 



X 
~X~ L 



4a tr (Tr/2) r' 2 



7TCT/i(7r/2) J 

<T tr (7r/2) 



cr tr (n/2) h, 



afii-K/2) H 
8 cr tr (7r/2) h 

m 



a fl {^/2y 3wa fl (7r/2) H 



cos 2 rdr 
at h rn < (HH) 1/2 . 



(58) 
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FIG. 5: The surface reactance X (in units of the reactance Xl in the absence of the DC current) vs. dimensionless amplitude 
h m /H of the AC signal at H = 500#. 



We notice that reactance (|5q ) is U£ r (7r/2)<7//(7r/2) ^> 1 times greater than that in the absence of the DC current. 
This is a direct consequence of the effect of amplification of electric signal at the film boundary which was treated in 
the previous section (see Eq. ji^)). The presence of strong DC current in the sample also causes linear behaviour of 
the reactance in the region of small amplitudes. As shown in Fig 5, the reactance decreases monotonically within the 
region between asymptotics (58) and (pTj). 



VI. CONCLUSION 



Nonlinear interaction of electromagnetic waves with a strong DC transport current in a thin metal film leads to 
unusual physical effects due to specific, typical only for metals, magnetodynamic mechanism of nonlinearity. These 
effects have been studied by analyzing the nonlinear response of the film, which carries the DC current, irradiated 
bilaterally by electromagnetic wave. The interaction of the wave with the current results in nonanalytical behaviour 
of the AC electric field on the sample surface which is characterized by appearance of sharp kinks. The increase of 
the current is accompanied by a rise of the amplitude of oscillations of the electric field at the surface of the sample. 
This, in turn, causes to the growth of the imaginary part of the surface impedance of the conductor. 

The results obtained in this treatment are valid under certain applicability conditions. Firstly, the AC electric field 
A£(x,t) must be small comparing to the potential electric field Eo(t). It follows from formulae (|2p|), (40), and (|4l| ) 
that the quantities £ and A£ m , Eq.(Q), are of the same order. Therefore, to ascertain the restrictions imposed by 
the condition £ -C Eo(t), we can use quantity A£ m in the latter condition. The quantity A£ m should be much less 
than the minimum value the function Eo(t), i.e. the magnitude of potential field ( |34| ) for cosujt = 0. The desired 
inequality reads 



2 h m l -2 / \ 
~HR n ^'' 



(59) 



where 5 n (u>) represents the characteristic penetration depth of the AC field into a metal under the condition of normal 
skin effect. Secondly, the non-uniform component of magnetic field inside the film must necessarily be much less than 
h m . This stems from the assumption that the AC magnetic field h(x, t) should be quasi-uniform (h(x, t) ~ h m coswi) 
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across the bulk of the film. The maximum value of the non-uniform correction can be estimated from the first 
of Maxwell's equations (||) as (ATra tr A£ m d/c) ~ h m (d/S) 2 , where an effective penetration depth 6(u>) equals to 
8 n (u})(R/l) 1 ^ 2 . As a result, we come to a requirement of the quasi-uniform property of the AC magnetic field which 
can be written in the following form: 

d 2 ^ (60) 

Comparing the restrictions imposed by inequalities (|59| ) and ( |6C| ) , it can be seen that condition (|5^) is more strict at 
large AC amplitudes, h m > H, while for small values of h rn one should use inequality (|60[). 

For a sample with thickness d — 10~ 3 cm, the electron free path I = 10 _1 cm, the concentration of electron N = 10 23 
cm' 3 , the Fermi momentum pp — 10 -19 g-cm/sec and for magnetic fields h m = H = 100 Oe, we have ui < 10 5 sec -1 
using conditions (|5^) and (|60|). At such values of the parameters, conditions J59| ) and (|6C]) are fulfilled as well as 
condition ([)]), which states that the value of the mean free path of electron should be large. 

Unusual manifestation of specific magnetodynamic mechanism of nonlinearity discussed in the present paper calls 
for future investigation. In particular it would be very interesting to explore experimentally theoretical predictions 
made in this work. 
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